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Abstract 

Usually, tc^tt^ decays are expressed in terms of weak amplitudes explicitly depen- 

dent on the CKM weak phase a or 7. In this letter, we show that the weak amplitudes 
can be rewritten such that a manifest dependence on (3 emerges instead. Based on this, 
we constrain new-physics contributions to the CP-violating phase (pa in B^-B^ mixing. 
Further, we apply reparametrization invariance and use QCD factorization predictions to 
investigate the bounds on an additional new-physics amplitude in B^ vr+vr". 
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One of the greatest successes of the B factories BaBar and Belle is the precise deter- 
mination of the CP- violating phase (pd in B mixing [Ij. In the Standard Model (SM), and 
using the Wolfenstein parametrization of the CKM matrix, 0^ is related to f3, one of the 
angles of the unitarity triangle, as (pd = 2/3. As B mixing is a loop process, the experi- 
mentally determined angle 0^ might in fact not equal 2j3, but be polluted by the effects of 
new-physics (NP) particles propagating in loops, thereby contributing an additional CP 
violating phase, see for instance Ref . [2] Q It is therefore of considerable interest to study 
any methods by which one can constrain an additional NP contribution to 0^. In this 
letter we shall show that the process B^ 7r"'"7r~ can be used to this effect. 

The set of neutral and charged B —>■ tttt decays has been extensively studied as a means 
of determining the angle a (or 7) of the unitarity triangle. The lack of a theoretically clean 
calculation of the strong amplitudes and phases involved can be overcome by exploiting 
isospin symmetry, see Ref. [Ij, commonly referred to as the Gronau-London method. It 
involves relating the various experimental observables (branching ratios and CP asym- 
metries) in i? — * TTTT to extract both the hadronic amplitudes determining these decays 
and the weak phase a. As an alternative to isospin, and in order to avoid B 
decays, the use of U-spin has been explored in Refs. j5] to extract 7 from 5° tt+tt^ 
and the U-spin related decay Bg K^K~ . In a conceptionally different approach the 
relevant strong amplitudes are calculated (as opposed to extracted from experiment), us- 
ing QCD factorization (QCDF) P [71 El E] or effective field theory methods (SCET) 
The advantage here is that less experimental input is needed, the disadvantage that the 
calculation is performed in a limit of QCD where the h quark is assumed to be very heavy. 
In any case, all these analyses put the emphasis on constraining the angles 7 or a. 

In this letter, we show that it is possible to express the decay amplitude in terms of 
(pd and /9, without any explicit reference to the angles a or 7. In the SM, the resulting 
expression allows the extraction of the relevant hadronic parameters from B^ -k^-k^ 
data alone, which can be compared to the theoretical calculation in QCDF. Beyond the 
SM, we study the possible presence of NP in this decay, which might contribute through 
mixing via a NP contribution to (pd 01 through an additional NP amplitude. 

We begin with a reminder of the parametrization used to extract 7. The amplitude 
for B^ — s> TT+TT^ is given in the form: 

^(5° ^TT + TT-) = A,A + A„A„, (1) 

where \q = V*^V^i,, and Ac, are strong amplitudes. A^ is dominated by tree diagrams, 
whereas the only contributions to A^ are from penguin diagrams. The corresponding 
time-dependent CP asymmetry is given by: 

^ r(g°(t)^7r+7r-)-r(g"(t)^7r+7r-) 

r{BO{t) ^71 + 71-) +r{BO{t) ^71 + 71-) 

= C± cos(Amt) — S± sin(Amt) . 
The experimental observables C± and S± can be expressed in terms of X±: 



2Im(A±) ^ 1-|A± 
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^NP at tree level this is highly disfavoured as its impact on B mixing would exceed the observed effects 
by far - unless the mass scales involved are in the ~ lOTeV range, see Ref. 
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Parameter 


Value 


Source 




X 


n 90^7+0.0009 
V.z,z,0 1 _o.0010 


PDG [UJ 








(41.2 ± 1.1) X 10-3 


PDG [m 




Vub 




(3.93 ±0.36) X 10-=^ 


PDG [m 




(21.1 ±0.9)° 


HFAG [1] 




(23.9 ±3.3)° 


this paper, Eq. (JT6i) 




7 


(77^i)° 


PDG [m 




Rb 


0.412 ±0.039 


this paper, Eq. 




Vul 

Vts 




0.214 ±0.005 


m 




Rt 


0.928 ±0.024 


this paper, Eq. (|T0|) 



Table 1: CKM parameters used in this letter. 



where X± is given by 



Ah 



A{BO -> vr+vr-) 
Parametrizing the amplitudes as in Eq. ([1]), we have 



re 



re' 



with re"^ = Ac/{AuRb), 7 = arg(-Ac/A„) and 



Rh 



A„ 



A. 



1-AV2 IK 



ub 



A 



IK 



cb\ 



(4) 



(5) 



(6) 



Numerical values for these and other CKM-related quantities are collected in Tab. [TJ The 
observables S± and C± are given by 



^4 



C4 



sin(0d ± 27) — 2r sin(0d ± 7) cos 5 ± sin 0^ 
1 — 2r cos 7 cos 5 + 

2r sin 7 sin 5 



(7) 



1 — 2r cos 7 cos 5 ± 



In the absence of penguin contributions, r = and the determination of (pd ± 27 would 
be completely analogous to that of (j)d from 5° J/ipKs- Realistically, r is expected to 
be a small, but non-zero number, which makes the extraction of 7 more challenging. 

We now show how a different parametrization of the decay amplitude ([1]) replaces 
the explicit dependence of A± on 7 by one on (3. Using [3 = arg(— At/Ac), one can trade 
the dependence on 7 for one on j3 by exploiting the unitarity of the CKM matrix and 
exchanging A„ for —Ac — A^: 



TT TT 



XcBc 
Ac(-Bc 



XtBt 

- Rte'^Bt) 



(9) 



2 



where Br 



-A,, and 





1 


\Vtd\ r 


Ac 


~ A 


\Vts\ I 



-(1 -2i?feCos7) A2 + 0(A^) 



(10) 



Note that B^ and Bt are both dominated by tree- level decays as they both contain Au- 
With this parametrization of the decay amplitude, A± becomes 



A 



± 



1 - RtRtce'f" 
1 - RtRtce-'P 

1 - de'^-^e'^ 
1 - de^de-i(^ 



(12) 



where Rtc = Bt/Bc and d = \RtRtc\, 0^ 
suppressed, but expected to be of order 1 
observables in ([3]) now read 



= aig{RtRtc)- Note that unlike r, d is not 
l^as Rt is also close to 1). The CP- violating 



^4 



C4 



^2 sin(2/3 



+ 2d cos 9d sin(0rf — /3) — sin( 



(i^ — 2d cos /3 cos 9d + ^ 
2d sin /5 sin 9d 



(13) 



(14) 



(i^ — 2d cos /5 cos 6^^ + 1 
Obviously ([T3|) . (|T^ are not independent of ([7]), ([8]), but related by the unitarity constraint 

Rte'^ + Rbe-'^ - 1 = . (15) 

The advantage of expressing S± and C± in terms of (3 instead of 7 is that, at least in the 
SM, there is now only one manifest weak phase. This implies that, with Rt determined 
from B mixing, both d and 6d can be extracted from experiment and compared to theo- 
retical calculations, for example QCDF. This is independent of any information from the 
decay B^ tt^tt^ whose branching ratio continues to be difficult to understand in the 
framework of QCDF or SCET. 

The most accurate measurement of (pd is via mixing in B^ decays to CP eigenstates 
of charmonium. The CP asymmetry averaged over these channels provides a direct mea- 
surement of sin^rf = 0.673 ± 0.023, so that in the SM (3 = (21.1 ± 0.9)° [10. It is also 
possible to derive f3 from tree-process measurements only, based on 7 and \Vub\- Taking 
7 and \Vub\ from Ref. [TT], see Tab. [H we find Ptree using 



sin Pi 



Rb sin 7 



tree 



2Rb cos 7 + Rl 



cos Pi 



1 — Rb cos 7 



tree 



^l-2RbCosj + Rl ' 



(16) 



which results in Pt^ee = (23.9^3;2)° (in the following analysis we use Ptree = (23.9 ± 3.3)°). 
Both values of P are compatible with each other, but we will use the latter one to obtain 
constraints on a NP contribution to (pd- 



^There is an ambiguity in this result, as /? = 
at the 95% confidence level by a Dalitz plot analysis of -B" 
time-dependent angular analysis of i?" ^ J/ipK*^ |14j . 



.9 ± 1.0)° is also a solution. However, this is excluded 
ZJ^/i^ where = 7r°, uj, 77 [13^, and by a 
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Experiment 


s± 


c± 


BaBar [15j 
Belle [IS] 


-0.68 ±0.10 ±0.03 
-0.61 ±0.10 ±0.04 


-0.25 ±0.08 ±0.02 
-0.55 ±0.08 ±0.05 


HFAG [IJ 


-0.65 ±0.07 


-0.38 ±0.06 



Table 2: Experimental results for S±, C± from BaBar and Belle and the HFAG average. 





I3k^,,s = (21.1 ± 0.9)° 


A,,, = (23.9 ±3.3)° 




d 




d 




BaBar 


0.790 ±0.031 


0.068 ±0.025 


0.775 ±0.037 


0.075 ±0.028 


Belle 


0.803 ±0.033 


0.158 ±0.041 


0.789 ± 0.040 


0.174 ±0.049 


HFAG 


0.796 ±0.021 


0.104 ±0.020 


0.782 ±0.027 


0.115 ±0.025 


QCDF 




-0 021+°-°''^ 


825+°-°^^ 

U.OZJ_Q Q52 


-0 021+°-°'^^ 

U.UZ1_Q Q42 



Table 3: Comparison of d and Od derived from experimental results to QCDF predictions for 
Pb^ccs from b — >■ CCS transitions and /3tree derived from tree decays. 

From the experimental data collected in Tab. [21 we find the values of d and 6^ given 
in Tab. [31 The high quality of the experimental results leads to small uncertainties on d, 
typically 5%, and moderate uncertainties on 6^. As expected, d is of order 1, while 6^ is 
rather smalljfl 

How do these results compare to theoretical predictions? As mentioned earlier, QCDF 
provides us with a framework for computing the individual amplitudes contributing to 
TT+TT^, and the relative phases between them. Complete results for both tree and 
penguin contributions are available to NLO accuracy [7J. The result for AJA^ isE 

4^ = -0.122^^1 - 0.024l°^^z. (17) 

Therefore using = A^ — A^, Bt = —A^, we find 

i?,. = | = 0.89n°^^-0.019l°S. (18) 

As Rt is, by definition, a positive number, 9d is given by the phase of Rtc, as seen in 
Tab. [3l The theoretical prediction is ca. 3a away from the experimental (HFAG) result 
and has the opposite sign. From Eq. (fTll) it is clear that 9d is mainly constrained by C±. 

■^Actually, there is a discrete ambiguity in the determination of d and 0d which yields a second solution 
d ^ 0.3 and 6d ^ 1- We discard this solution as d is the ratio of two tree-dominated amplitudes and hence 
expected to be close to 1; d i-^^ 0.3 would imply a massive NP amplitude which can not be generated in 
any NP models we are aware of. 

Recently the NNLO calculation of the tree amplitude has been completed in Ref. [S] . We do not use 
this result in our analysis, as we also require the penguin amplitude at NNLO. Nevertheless, we have 
checked numerically that using this value would alter Ac/A^ by 4 to 5%. 
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While Tab. [2] clearly shows that the results from BaBar and Belle are not completely in 
agreement, they do agree on the sign and require a positive 6d- This situation is similar 
to that for B Kir decays, where the sign of the observed direct CP violation is difficult 
to reconcile with QCDF predictions. The source of the discrepancy is independent of 
any NP contributions to Arrid^s, from which Rt is determined, as 6d is independent of 
Rf. It is also unrelated to any NP contributions to the B mixing phase 0^, as using the 
tree- level determination of /?, Ptree, for the explicit P dependence renders the discrepancy 
even slightly worse (right-hand side of Tab. [3]). An explanation in terms of NP could 
be that an additional NP amplitude contributes to B^ —>■ tt'^'tt", and we explore this 
possibility below. Otherwise, the discrepancy may be related to neglected higher order 
terms in QCDF - either from radiative corrections or terms suppressed by inverse powers 
of the b quark mass. As for extreme values of the QCDF input parameters a positive 6^, is 
possible, it would be interesting to see what choice of these parameters yields positive 
and whether this choice agrees with the scenarios advocated in Refs. [3 El E] to reconcile 
the QCDF predictions for B{B^ — > tt^tt^) with experimental data. 

In order to calculate d from QCDF we also need Rt- This can be determined from 
the ratio of the mass differences in the Bg and B^, systems. CDF has made a very clean 
measurement of ArUg [17] , which can be turned into a value for Rt using lattice information 
on the relevant hadronic parameters [12j and Eq. (fTOl) . The result is given in Tab. [H This 
value of Rt allows for a limited contribution of NP to the mass differences in the B^^s 
systems: from the unitarity relation (ITSi) . using Rt and (3 as input, one finds Rb ~ 0.36 
and 7 ~ 68°, in good agreement with the values quoted in Tab. Combining the QCDF 
calculation with this value for Rt, we find that the QCDF prediction agrees with the 
experimental results within la. The experimental results and QCDF predictions for d 
and 6d are shown in Fig. [TJ The analogous QCDF predictions needed in the "standard" 
parametrization of 5*^, Eq. ([7]), and C±, Eq. (jS]), can be found by dividing (ITTj) by Rh- 

r = 0.30lJ^, 5 + 7r = 0.19m^ (19) 

Note the marked difference in the relative errors for {d, 6^) and (r, 6), which is an impor- 
tant advantage to the parametrization (Q of the B^ vr+vr" amplitude and will reduce 
the uncertainties in the constraints on a NP amplitude contribution to this decay. 

Let us now discuss constraints on 0^ — 2/?, based on the QCDF results in Tab. [31 
Assuming there is no NP amplitude contributing to the decay, one can obtain f3 using 
(pd and d, 6^ from QCDF as input. This results in f3 = 19.9° with a minimum of 
63. Leaving 6^ as fit parameter instead, we find (3 = (16.4 ± 6.5)° which shows that 
the above discrepancy in 6d from QCDF and from experiment is largely irrelevant for 
NP contributions to (pd- Both values are consistent with 0^/2 = (21.1 ± 0.9)°, and and 
also with (3 = (23.9 ± 3.3)° determined from 7 and \Vub\- However, it is still interesting 
to note the two results for 0^/2 — j3 are of opposite sign. If there was indeed a NP 
contribution to 0^, we would expect to see common trend in the sign of 0d/2 — j3 from 
various determination. As this does not seem to be the case, we take this as an indication 
that any NP contribution to 0^ is indeed small. 

Given this result, the decay B^ tt'^tt~ seems then a very suitable place to constrain 
any NP contributions to the decay amplitude. Such contributions are possible in a variety 

^The agreement is not perfect for If, however, one uses the result for \Vub\ ~ 3.5 x 10^"^ from 
exlusive decays, instead of the larger PDG value, Rb ~ 0.36 and the agreement is near perfect. 



5 



0.25 
0.20 
0.15 
0.10 
0.05 
0.00 
-0.05 




0.75 0.80 0.85 
d 



0.90 



0.25 
0.20 
0.15 
0.10 
0.05 
0.00 
-0.05 




0.75 0.80 0.85 
d 



0.90 



Figure 1: The allowed experimental region for d and 9d are plotted for the case f3 = (21.1ib0.9)° 
on the left and (3 = (23.9 it 3.3)° on the right. The blue (lower) region corresponds to the results 
for S± and C± of BaBar, the green (upper) region corresponds to those of Belle and the red 
(central) region corresponds to the average given by HFAG as in Tab. [2j The blue point with 
error bars represents the QCDF prediction. 



of NP models, e.g. the 2-Higgs doublet model [TS] or R-parity violating SUSY [T^. We 
therefore introduce the (complex) NP amplitude Anp and a new weak phase 5np, such 
that 

\1 - de^de-il^ + A^pe-^^^p J 

Here d and 6^ are SM quantities, given by the QCDF values stated earlier, in Tab. [31 The 
above expression for A± is actually not suitable to constrain A^p and (5np: in Ref. [20], it 
was shown that a given amplitude, with an arbitrary number of distinct weak phases, can 
always be expressed in terms of any two weak phases. Using this so-called reparametriza- 
tion invariance, A± can be expressed in terms of the two weak phases (f)d and 

where 

^,^ie'^ ^ o^e^^'^ - ANpsin((5Np)/ sin/3 
1 + /Inp sin(/3 - 5np)/ sin /? ■ 

Note that d'e'^^'^i d'e^^'d under a CP transformation, so 9'^ is indeed a strong phase. 
In Eq. (12T|) . 0^ is, by definition, the phase measured in B mixing, i.e. from h — > ccs 
transitions, while j3 is obtained from tree-level processes, i.e. ipd = (42.2 ± 1.8)° and 
f3 = (23.9 ± 3.3)°. The values for d' and 6'^ are obviously identical to the experimental 
results for d and 6d in Tab. [31 The resulting constraints on Anp and 5np are shown in 
Fig. [21 Depending on the value of (5np, large NP contributions |Anp| are possible. Note 
that these results are mainly due to the discrepancy between 6'^ from experiment and 6d 
from QCDF. In particular, even for 5np = a non-zero |v4np| ~ 0.1 is needed. Also note 



^Eq. (HH) was already obtained in Ref. [10], as (48). 
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Figure 2: |^np| and arg(^Np) plotted against the weak phase of the new physics amphtude, 
(5np- The thick black lines represent the results found using the experimental average value of 
sin(^ci), the grey regions represent the uncertainty. 



that the allowed valued for |Anp| can be as large, or even larger, than the SM (QCDF) 
prediction for d. Although it would be interesting to interpret the NP amplitude in terms 
of the NP models mentioned, making use of constraints from other processes, this is 
beyond the scope of this letter. 

In summary, we have investigated the effect of NP on the decay -B° — > tt+tt". We first 
analyzed the CP asymmetries in the SM, using a particularly convenient parametrization 
of the decay amplitude which depends on only one weak phase, (3. We compared the 
experimental results of the relevant hadronic quantities with those from theoretical calcu- 
lations, using QCDF. We found that while the size of the strong amplitude is consistent, 
its predicted phase 9d deviates from experiment by ~ 3a. It would be interesting to study 
whether the choice of QCDF input parameters minimizing this discrepancy is the same 
needed to reduce that between the predicted and observed branching ratio of S° — > 7r°7r°. 
Wc then analyzed the implications for NP contributions to the decay. Under the assump- 
tion that NP modifies only the phase in B mixing, the discrepancy in the strong phase 
was largely irrelevant and the NP contributions to the mixing phase found to be consis- 
tent with 0. We then allowed for a NP amplitude contributing to the decay amplitude 
of S° — > 7r+7r~ and analyzed the CP asymmetries exploiting the constraints imposed by 
reparametrization invariance. Here we found scope for large NP contributions, which are 
driven by the discrepancy between 9^ from experiment and from QCDF. An analysis of 
this NP amplitude within specific models is beyond the scope of this letter. 
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